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Abstract 

We analyze the dynamics of a nanomechanical oscillator coupled to an 
electrical tunnel junction with an arbitrary voltage applied to the junction and 
arbitrary temperature of electrons in leads. We obtain the explicit expressions 
for the fluctuations of oscillator position, its damping/ decoherence rate, and 
the current through the structure. It is shown that quantum heating of the 
oscillator results in nonlinearity of the current- voltage characteristics. The 
effects of mechanical vacuum fluctuations are also discussed. 
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I. INTRODUCTION 



The rapid development of nanotechnology in recent years has ushered in a new generation 
of quantum electronic devices incorporating the mechanical degrees of freedom [|l]||, so- 
called nanoelectromechanical systems (NEMS). A particular example of such a device is the 
tunnel junction having its transition matrix element modulated by a vibrational motion. 
This modulation can be achieved by the introduction of small grains (shuttles) embedded 
in the elastic medium between the leads [0,(|,|5],|6],|7|] , or by coupling the tunnel junction to 
a mechanical oscillator (cantilever) flj9|,p!0|]. These two cases can be described by similar 
formulations but the oscillator masses are very different (for example, m ~ 10 -20 g for the 
C*6o molecule used as a shuttle and m ~ 10 _12 g for the cantilever). Furthermore, in addition 
to nanomechanical applications, such systems have a direct relation to the fundamental 
quantum theory of measurements fLl] , |12| . In this context, the tunnel junction can be consider 



as a device measuring the position of the oscillator. As was shown in Ref. |10| for the case 
of a large voltage applied to the junction, this measurement process induces decoherence 
and dephasing of the oscillator, even at zero temperature. From another point of view, the 
nanomechanical oscillator can serve as a measurement device for the electronic subsystem. A 
similar situation takes place in magnetic resonance force microscopy |13| , where the magnetic 
moment of electron or nuclear spins in the sample (attached to the cantilever) is measured 
by means of optical interferometry of the cantilever motion. 

There are two manifestations of the effect of mechanical oscillations on system behavior 
and, consequently, there are two characteristic scales. The first phenomenology occurs when 
the uncertainty of the oscillator position is of the order of the tunneling length. This is 
always the case for the shuttle system, but for the cantilever it only exists at high bias 
accompanied by oscillator heating. (In the present paper we simplify the model used in 
Refs. ||[|||, such that our shuttle system has a permanent electrical connection to one lead, 
while the connection to the other lead is mediated by the tunnel junction.) 

The second class of phenomena is associated with resonance between the oscillator fre- 
quency and the bias voltage applied to the tunnel junction. The characteristic frequencies of 
the oscillator are of order of 1GHz |T|JT^1 and this resonance occurs at sufficiently low bias. 



It should be noted that the theoretical approach used by Mozyrsky and Martin in Ref. [10 
fails at low bias voltage and requires improvment. 

NEMS can be used as an ultrasensitive device for magnetic and biological applications 
and, therefore, its own fluctuation level - determining the signal/noise ratio - is of crucial 
importance. Moreover, quantum measurement procedures also require a high level of sensi- 
tivity. The mechanical oscillator can be considered as an open quantum system interacting 
with two electron reservoirs (leads) serving as an effective heat bath. In the present pa- 
per we employ the general theory of open quantum systems developed in Refs. p~5| , p76|JT7[] 
to determine the nonequilibrium fluctuations of the mechanical oscillator for arbitrary bias 
voltage and temperature. Current noise created by the junction is assumed to be the main 
mechanism of decoherence of the oscillator. 

We derive general formulas and employ them in the two cases described above, showing 
that at low temperature the bias applied to the tunnel junction leads to oscillator heating 
at voltages larger than the critical value associated with the characteristic frequency of the 
oscillator. With increasing temperature this resonant picture is smoothed and, moreover, 
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for higher temperature the fluctuation level approaches that of an uncoupled oscillator. 
Furthermore, we show that in the shuttle nonlinearity of the oscillator-junction 

interaction makes a pronounced contribution to the level of mechanical vacuum fluctuations 
even for zero bias and weak coupling between the above-mentioned subsystems. 



II. GENERAL FORMALISM 

The Hamiltonian of the system tunnel junction and mechanical oscillator is given by 

H = H L + H R + H tun + H . (1) 

Here, 

H a — X] Ek a ct a c ka (2) 
k 

is the Hamiltonian of the left, right leads (a = L, R, respectively) 

p 2 mujx 2 

is the Hamiltonian of the mechanical oscillator with a mass m and the resonant frequency 
u>o,, and 

H tun = - J2(T kq 4 L c qR + h.c.)e~ x ' x (4) 

kq 

is a tunneling term depending on a position of the oscillator, where A is the characteristic 
tunneling length and h.c. is the Hermitian conjugate. The electron gas in the leads plays 
the role of a heat bath for the oscillator and their nonlinear interaction given in Eq.(4) can 
be written in the form 

Htun = -Qe-*' x (5) 
where Q(t) is interpreted as an effective heat bath variable given by 

Q{t) = Y^{T kq c + kL c qR + h. C .). (6) 

kq 

The heat bath thus defined is characterized by a response function ip(t, t±) and a symmetrized 
correlation function M(t,ti) of the unperturbed variables Q^(t) taken in the absence of 
interaction (h = 1, k B = 1), as 

¥> (*,t 1 ) = <<[Q (0) (*),Q (0) (*i)]->«(*-*i), 

M(t,t 1 ) = ^[Q (0) (*),Q (0) (*i)] + >, (7) 

where 6{r) is the unit Heaviside step function and [...,...]+ and [...,...]_ are the anticom- 
mutator and the commutator, respectively. The angle brackets refer to averaging over the 
equilibrium states of both the left and right leads at the same information. The chemical 



potentials of the leads can be different with Hr — — eV, where V is a voltage applied to 
the tunnel junction. For weak coupling between the mechanical oscillator and the electronic 
bath (weak tunneling) the action of the oscillator on the dissipative environment is described 
by the formula [15] 



Q(t) =QW(t) + J dtMt,h)e 



x(ti)/A 



(8) 



In turn, the effect of the dissipative environment on the mechanical oscillator is determined 
by substituting Eq.(8) into the Heisenberg equation of motion for the position operator of 
the oscillator given by 



X + LU X 



Q(t)e 



-x/X 



(9) 



To eliminate the unperturbed heat bath variables Q^°\t) from Eq.(9), we apply the quantum 
Furutsu-Novikov theorem [|T5] 



(Q®(t)e-*W X ) = f dti(g^(t),g (0) (*i)>< 



-x(t)/X 



(10) 



where 8/8Q^°'(ti) is the functional derivative with respect to the free heat bath variable 
Q(°>(ti). The functional derivative of an arbitrary operator A(t) of the dynamical system is 
proportional to the commutator |15| in the form 



6A(t) 



% r 

h 



A{t),e 



-x(ti)/X 



9(t - h 



(11) 



As the result, we obtain the non-Markovian stochastic equation of motion for the position 
operator of the oscillator as given by 



X + UqX = f (t) 



Am 



dti M(t,ti)i 



e -x(t)/X e -s(ti)/A 



e *(*)A e -*(*i)A 



;i2) 



where M(t,t\) = M(t,t\)6(t — t\) is the causal correlation function of the heat bath and 
is the fluctuation source, 



g(o) (t)e -W/A_ f dt x {QW{t),Q®{h))i 



e -x(i)/X e ~x(h)/X 



eit-h)), (13) 



having zero mean value, (£) = 0, according to Eq.(10). The Langevin-like equation, Eq.(12), 
as well as the whole method of Refs. []15,|16,r7|, takes into account a nonlinearity of coupling 
between the subsystems and, in addition, incorporates the nonlocal character of heat bath 
fluctuations. In these respect this treatment goes beyond the well-known Caldeira-Leggett 
approach [ p~8|| . 

In the case of weak coupling, the correlation function of fluctuation sources is given by 



<2 [£(*U(*')] + > 



^^{({Q w (tlQ (0) (t')} + )( 
+([Q (0) W,Q (0) (O]-)( 



e -x(t)/X e -x(t')/X 



) 



e -x(t)/X e -x(t')/X 



)}• 



(14) 



A 



It should be mentioned that the fluctuations of electronic variables {Q^(t)} are non- 
Gaussian. However, as indicated above, Eqs.(8),(10),(14) are valid for weak coupling between 
the dynamical system (mechanical oscillator) and the dissipative environment (electrons in 
the leads). In this case we can calculate the (anti) commutators in Eqs.(12),(14) using a free 
evolution approximation 



x(t) = x(ti) cosuj (t — ti) + 1—L- sinuj (t — ti), 

TTLUOq 



and employing the Baker-Hausdorff theorem, we obtain 



e -x(t)/X e -x{ti)/X 



e -x(Jt)/\ -x(Jti)/\ 



cos(f sintuor) exp 



2 sin(z/ sin uqt) exp 



x(t) +x(ti) N 

x(t) + x(t 1 y 

A 



(15) 



(16) 



Here, 



h 



I/ = 



2muJoX 2 



(17) 



is the square of the ratio of the unperturbed oscillator position uncertainty at zero tempera- 
ture and the tunneling length A. The operator x(t) can be written as a sum of the mean and 
fluctuating parts, x(t) = x(t) + x(t). Even if the averaged position of the oscillator is small 
compared to the tunneling length, x(t) <C A, the fluctuating amplitude can be of the order 
of the tunneling length causing nonlinear effects. To examine this nonlinearity, we assume 
that oscillator position fluctuations are approximately described by Gaussian statistics with 
a dispersion (x 2 ) = (x(t)x(t)). In this approximation, we obtain the exponent involved in 
Eq.(16) as 



exp 



xjtf + xjtiY 



1 - 



A 



A 



exp 




+ ((l/2)[x(Q, a-fo )]. 
A 2 



Eqs.(16) and (18) provide the expressions required on the right side of the non-Markovian 
stochastic equation, Eq.(12). The dispersion (x 2 ) and the correlator ((l/2)[x(t),x(ti)]+) 
have to be determined self-consistently. To accomplish this, we suppose that the operator x 
of the oscillator position obeys the equation 



X + "fX + CUqX 



(19) 



where the damping rate 7 can be calculated from the collision term on the right of Eq.(12). 
Using Eq.(19), we obtain correlator and the dispersion of the oscillator position as 



<5 2 > 



K(u ) 



2^o7 ' 

(5 2 ) e - 7l *-* l|/2 coawoit-h). 



(20) 



Here, 

K{u>) = Jdre^{\Ur),m\ + ) (21) 

is the spectral function of the fluctuation forces, Eq.(14). Assuming that the averaged 
deviation of the oscillator position is much less than the tunneling length A, i <C A, (the 
root-mean-square amplitude of mechanical oscillations can be of the order of A), we obtain 
a simplified equation for the total coordinate x = x + x as 

x + uolx + I dhL(t - h)[x(t) + x(h)} = £(t) - —F , (22) 

where the collision kernel L(r) is given by the expression 

L(t) — 2M(r) sin(z/ sinc^oT) + v?(r) cos(z/ sincij r) exp[z/ c (l + coscj t)] (23) 
with the parameter 

"c = ^ (24) 

describing nonequilibrium fluctuations of the oscillator position. The deterministic force F 
on the right side of Eq.(22), which arises from interaction of the oscillator with the tunnel 
junction, is given by F = (l/A)A(u; = 0), where A (a;) is the Fourier transform of the kernel 
L(t), Eq.(23). It follows from Eq.(22), that coupling to the tunnel junction also results in 
a shift of the oscillator frequency, cJq — > — (l/mA 2 )A(cj = 0), jointly with damping and 
decoherence described by ImA(ui) = A" (a;). 

The relaxation rate of mechanical oscillations due to coupling to the tunnel junction is 
given by 

7 = 4^*1, (25, 

where 

1 i = + oo 
l=— oo 

j i=+oo n=oo 

+ o e ^ E £ M^^nKHx'V + 2n^ - Zw ) + x'V - 2nw + lu ) 

l=—oo n=l 

+x"(^ + 2nw + Ivq) + x"{u ~ 2nuj - Ilu )} 

i i=+oo n=oo 

+2 e " E E W-WnO 

£=— oo n=0 

x{S'[cj - (2n + l)u; - ^o] + ^[a; - (2n + l)u + Zcj ] 
-Sfw + (2n + l)w - luo] ~ S[lo + (2n + l)u + Zw ]}- (26) 

Here, J n (v) is the Bessel function of n-th order, and hiy) is the modified Bessel function of 
1-th order. The spectral function S(u>) and the imaginary part of the susceptibility x"^) are 



the Fourier transforms of the correlator and the response function of the heat bath, Eq.(7), 
respectively, 

S(u) = J dre iu}T M(r), 
X» = JdrsmicorMr). (27) 

For a harmonic oscillator, the inverse relaxation rate 7 _1 given by Eq.(25) also represents a 
characteristic time scale for the coherence decay, Td = 7 _1 . 

The parameter v c of Eq.(24) describing the level of mechanical fluctuations can be found 
from the self-consistent equation 



KM 
Vc 2cu 2 AV 

where K(u) is given by 

*2 Z=+oo 

%) = ^ C E ii(yc)Jo(yo){S(w-M + S{w + M} 

l=—co 

j i=+oo n=oo 

+ -e^ E E / i( z/ c)^2n(t / o){5'(^ + 2ncu - lu ) + ^(a; - 2n^ + M 

l=— 00 n=l 

+S(u> + 2nuo + Iujq) + S{u — 2nu — lu>o)} 

■y l=+oo n=oo 

-o e " c E E JiOO-WiMfa'V + (2n + 1 Vo - Zw ] - X'V - (2n + l)w - Zw ] 

2 /=-oo n=0 



(28) 



+x'V + (2n + l)w + /cj ] - - (2n + l)w + ^o]}- (29) 



III. RESERVOIR CORRELATION FUNCTIONS AND TUNNEL CURRENT 

To determine the oscillator damping rate 7 of Eq.(25) and the level of fluctuations v c of 
Eq.(24), we need to determine the spectral function S[iS) and the susceptibility }("{oj) of the 
electronic heat bath, Eq.(27), for the case with no interaction with the oscillator. Assuming 
that the variables {cfl, c^, ..} of the free electron gases in the leads obey the Wick theorem 
(with uncorrelated electron systems in the different leads), we can express the correlator of 
unperturbed heat bath variables in the form 

<Q<°>(f)Q (0) (*i)> = E \T kq \ 2 {{cf L + (*)cS(t 1 )>( C W(t) c W+ (tl)> 

kq 

+ (cS(t)cS + (ti))(4°i + w4°i(^))}- (30) 

Introducing retarded, advanced and " lesser" Green's functions of electrons in the leads as 

sL(Mi) = ((-<)[cS(0,ci2 + («i)] + >^-«i) = -ie- iE ^-^9(t- tl ), 

<?L(Mi) = (<[cS(0,cS2 + (*i)] + >f(*i ~t) = ie-^-'Mh-t), 

ff£(Mi) =i(c$(t),cV> + (t 1 ))6(t-t 1 ) = if a (E ka )e- iE ^\ (31) 
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where 



f a (E) = f(E-fi a ) 



. E - fj, a 
exp — — — + 1 



-l 



(32) 



is the Fermi distribution in the a— lead (a = L, R) with the chemical potential \i a and 
temperature T, we obtain 

M(t,h) = l^lT^g^hM^t) -g a kL (ti,t)} 

kq 

+[^fl(Mi) - # fl (Mi)te^ 

+[9 r kL (t,ti) - <72 L (Mi)]«7?*(*i,*) + 2g< R {tM)gt L {tut) + 2g< L {tM)gU t ^ t )}^ 

<p(t, h) = i6{t - t ± ) \T kq \ 2 {[9 r gR (t, h) - g a qR (t, hM^t,, t) - g< R (t, hM^t) - g a kL {t u t)\ 

kq 

+[9 r kL (t,ti) - g a kL (t:ti)}g< R (t 1: t) - g^ L {tM)[g r qR {ti,t) -g a qR (t u t)]}. (33) 

As a result, the formulas for the spectrum and the imaginary part of the susceptibility of 
the dissipative environment are given by 

S(u) = tvJ2 \T kq \ 2 [S(uj - E kL + E qR ) + 5{u + E kL - E qR )\ 

kq 

x [f R {E qR ) + f L (E kL ) - 2f L (E kL )f R (E qR )}, 
X » = tt]T \T k q\ 2 [S(uj + E kL - E qR ) - 5(cu - E kL + E qR )][f L (E kL ) - f R (E qR )}. (34) 

kq 

It should be mentioned that when the leads are in the full thermodynamic equilibrium state, 
(fii = Hr), the spectral function S(u>) and the function are related by means of the 

fluctuation-dissipation theorem S(cu) = x"{oS) coth(cj/2T). 

The electric current through the junction is defined as I — II — —Ir, where I a = e(N a ) 
with N a = J2k c ta c ka being an electron number in the a— lead. It follows from the equations 
of motion for the electron operators of the leads, 

tCpL = E pL c P L — ^2T pq 'C q ' R e~ x ^ x , 

q' 

ic qR = E qR c qR - T p ' q c p > L e- x/x , (35) 
k> 

that the electric current depends on the oscillator position x as 

I = -*£<{^c 9 VpL - h.c.}e- x ' x ). (36) 

pq 

In the case of weak tunneling between leads, the electron operators of the leads can be 
represented as 

c P L(t) = cf L {t) -]TT P(? , J dhg^h^i^e-^, 
q' 

c qR {t) = 4°i(0 - ET p %/ *i^,«i)cJ?l(*i)e- x(tl)/A . (37) 



8 



Therefore, the electric current is given by 

I = E \ T v q \ 2 I dti{g r pL (t,ti)g< R (t l7 t) 

pq J 

-9 r qR {tM)9^t)}{\ [e-°W\e-°W' x ] + ) + h.c. (38) 
Substituting the expressions for the anticommutator, Eq.(16), we obtain 

i=+oo n=+oo 

/= o e " E E ^K)J2n(^0)El^| 2 [/^)-^(^)] 
i=-con=-oo pg 

x {5[£ pL - E qR + 2n)u ] + 5[£ pL - E qR -{I- 2n)u ] 
+5[E pL - E qR + + 2n)iu ] + 5[E pL - E qR - (I + 2n)iu ]}. (39) 

We can determine the spectral function S(cu) and the imaginary part of the susceptibility 
x"(w) of the electron reservoir, Eq.(34), by introducing the densities-of-states D a (E) of the 
leads and replacing sums over k, q, .. by integrations over the corresponding energies as 



£(...) - J dE L dE R D L {E L )D R {E R ){...). 



kq 

We assume that the tunneling elements \T kq \ 2 do not depend on energies, \T kq \ 2 = |T | 2 , 
and the densities-of-states near the Fermi surface in the leads are also energy-independent, 
D a {E) ~ D a (/i) = D a with /i R — /i + eV/2,fj,L — — eV/2,a = L,R. Furthermore, we 
assume that all energy parameters of our problem (eV, T, uq, ..) are much less than the basic 
chemical potential \x of the electron gas in the leads. In this case, we obtain the imaginary 
part of the susceptibility as 

/+oo 
dE[f(E-uj)-f(E + uj)}, (40) 
-oo 

whereas the spectral function S(u>) has the form 

/+oo 
dE{[l - f(E - n L )][f(E -fi R + uj) + f(E -HR- uj)} + 
-oo 

f(E - n L ) [2 - f(E -fi R + uj)- f(E -HR- co)]}. (41) 

As a result, the electrons in the leads represent an Ohmic heat bath with respect to the 
mechanical oscillations of the cantilever with the imaginary part of susceptibility given by 

X» = au, (42) 
where a = 2n D L D R \T \ 2 , and the frequency- dependent spectral function has the form 

'uj — eV\ , Tr . , fuj + eV*' 



m = f 



i (0J-eV\ . Trx , fu + eV\ 
(u - eV) coth [-^T^) + + eV) coth (^^J 

At zero temperature S(u) is given by 



(43) 



S(u) = a[u9(u-eV)+eV9(eV-u)], (44) 

where 9{uj) is the Heaviside step function. 
For the electric current, we finally obtain 

/ = G(V)V = ae 2u °eV, (45) 

where G(V) = eaexp(2z/ c ) is the nonlinear conductance of the tunnel junction, which de- 
pends on the fluctuation level of the mechanical oscillator. 



IV. RESUTS AND DISCUSSION 



Although the general formalism developed in the previous sections makes it possible to 
perform comprehensive analyses of the rich phenomenology occurring in a wide variety of 
nanoelectromechanical systems, we have restricted ourselves in the present paper to the 
examination of nonequilibrium fluctuations of a mechanical oscillator coupled to a biased 
tunnel junction, along with the study of quantum heating effects on the decoherence rate 
(Eq.(25)) of the oscillator and on the current- voltage characteristics (Eq.(45)) of the junction. 

To start, we demonstrate analytical results related to the case of weak nonlinearity of the 
cantilever-junction coupling, Eq.(5), when u <C 1. In this case, the contribution of vacuum 
fluctuations to the collision term A" (a;), Eq.(26), and to the correlator K(u), Eq.(29), can 
be neglected. Therefore, the dissipative kernel is proportional to u, A"(c<j) = aexp(2i/ c )co>. 
This corresponds to Eq.(19) for the oscillator position with the damping rate 7, Eq.(25), 
and the spectrum of the fluctuation forces K(u), Eq.(29), given by 



n 



7 = a 



mX 2 



ti 2 



l=+oo 



2m 2 X 2 



e Uc E Ii(v c )[S(u-Iu )+S(u + Uj )]. 



(46) 



i=— 00 



In contrast to the approach of Ref. [10|, the fluctuations of the random forces £ involved in 



the Langevin equation, Eq.(19), are not white noise because of the frequency dispersion of 
the correlator K(u), Eq.(46). Even for the case of weak heating, when the root-mean-square 
amplitude of the cantilever fluctuations is much less than the tunneling length A, J (x 2 ) <C A, 
and u c <C 1, the spectrum K(u) of the fluctuation sources, 



K{u>) 



h 2 



a 



2m 2 X 2 



eV) coth 



uj - eV\ . . . fuj + eV 

+ (00 + eV) coth 



2T 



2T 



(47) 



demonstrates a non-white character, especially at zero temperatures of the electron reser- 
voirs. 

On the basis of the above-mentioned approximations and with Eqs. (28), (46), (47) we 
have derived the dispersion of the cantilever fluctuations (x 2 ) = X 2 u c as a function of voltage 
V applied to the junction and of temperature T of the electron reservoirs, as 



x 2 ) 



h 



4muj 



'u 



w V 2T . 



uj V 2T 



(48) 



1 n 



At low temperatures, T <C u ± eV, and at low voltages, eV < t^o, when the current in the 
tunnel junction can not stimulate an excitation of the mechanical system, the dispersion 
of the oscillator fluctuations remains on the vacuum level, (a; 2 ) = (h/2muJo). However, at 
the higher voltage, eV > u Q , the level of fluctuations increases linearly with the voltage, as 
(x 2 ) = (h/2mu> )(eV/u) ), and we find the dimensionless parameter v c as u c = u (eV/u ). 
This linear increase of the fluctuation level may be interpreted as an increase of an effective 
temperature of the oscillator: T e ff = eV/2, where T e ff is proportional to the dispersion of 
mechanical fluctuations, (x 2 ) = T e ff/moj 2 . This confirms the result obtained by Mozyrsky 



and Martin for the Caldeira-Leggett model. 

The heating process affects also the bias dependence of the damping rate 7, Eq.(46), 
and the behavior of the conductance G(V), Eq.(45). At low voltages, eV < Hlj , both of 
these characteristics do not depend on V, [G = ea(l + 2z/ ),7 = {ha/mX 2 )(l + 2h> )], whereas 
above the threshold of the excitation of cantilever oscillations, eV > huo, the heating process 
makes a linear contribution to conductance and to the relaxation rate: 

G(V) =ea(l + 2u ^-) , 7 = (l + 2u ^-) . (49) 



hu! J mX 2 \ hxii 

We consider here two models of an oscillator having different masses. For the first model, 
the shuttle, the mass is 10 _20 g, and this model describes the molecules embedded in an elastic 
medium between the leads. The molecule is taken to be electrically connected to one lead, 
whereas a connection to the other lead is realized by a tunnel junction. 

The second model is associated with the cantilever having a mass 10~ 12 g which is placed 
in the neighborhood of the electrical contact; in this case variations of the cantilever position 
modulate the tunneling matrix element of the junction. The oscillator frequencies are the 
same for both cases and are equal to 1GHz, hu = 10~ 18 erg. The tunneling lengths are also 
the same both for the shuttle and the cantilever and are equal to 10 _8 cm. As a result, the 
parameter u is much less than one (~ 5 ■ 10~ 7 ) for the cantilever and is of the order of one 
(~ 0.5) for the shuttle. 

To characterize the nonequilibrium fluctuations of the systems, we plot the bias and 
temperature dependencies of the parameter 

v c 2muj {x 2 ) 



coth(huo /2T) h coth(huo / 2T) 

i.e. the ratio of the dispersion fluctuations of the oscillator coupled to the tunnel junction 
to that of the uncoupled oscillator at temperature T. 

The bias dependencies at various temperatures are presented in Fig. 1(a) for the cantilever 
and in Fig. 1(b) for the shuttle. One can see that at extremely low temperature there is 
the change of the curve shape at the critical voltage corresponding to the characteristic 
frequency of oscillator. For the cantilever, there is not any heating below this bias and 
there is linear growth of fluctuations thereafter, as indicated by the analytical results above. 
As temperature increases the curve becomes smoother and heating decreases, and, finally, 
there is almost no heating at T = 0.1K. For the shuttle case, there is some heating even 
below the critical voltage and, moreover, the equilibrium value (at V = 0) of v c is not one. 
The reason for this is that at ~ 1, vacuum fluctuations in the leads (back and forth 
tunneling events) give rise to oscillator anharmonicity and the steady-state fluctuation level 

1 1 



of the oscillator nonlinearly coupled to the electron reservoirs is different from that of the 
uncoupled oscillator. 

The temperature dependencies are exhibited in Figs. 2, 3(a,b) and 4(a,b) for the can- 
tilever, and similar features have been found for the shuttle. The bias voltages are chosen 
to be below (Fig.2), above (Fig.3(a,b)) and far above (Fig.4(a,b)) that of the excitation 
threshold huj^/e = 0.625 x lO -6 ^, respectively. The "b" parts of the figures magnify the low 
temperature behavior. It is evident from these figures that as temperature increases, v c ap- 
proaches one for any bias voltage. However, for low bias there is a peak at low temperature 
which becomes more pronounced below the critical bias voltage. This peak occurs when 
temperature is such that T ~ u>o — eV; the peak is absent at high voltage. This phenomenon 
occurs for both the cantilever and the shuttle. 



V. SUMMARY 



In summary, we have applied the theory of open quantum system [^,^,|T^] to the na- 



noelectromechanical system consisting of an oscillator coupled to a tunnel junction. We 
have obtained explicit expressions for the nonequilibrium fluctuation level of the oscillator 
(Eqs.(28),(48)) as well as for the tunnel current through the structure (Eq.(45)) and for the 
decoherence rate of the oscillator (Eqs.(25),(46)). The bias and temperature dependencies of 
the oscillator fluctuation level have been determined. Considering two specific models with 
different oscillator masses, we have shown that, for small mass, the level of mechanical vac- 
uum fluctuations in this system differs significantly from that of a linear harmonic oscillator 
with the same frequency ujq because of an anharmonicity induced by the weak nonlinear 
oscillator-current coupling. This effect takes place even at zero bias voltage, but vanishes at 
higher temperatures. We have also found that for voltages below the excitation threshold, 
huoo/e, the relative level of oscillator fluctuations (normalized to the equilibrium value) peaks 
at temperatures of the order of the difference between the characteristic frequency of the 
oscillator and the applied bias. 
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FIGURES 



FIG. 1. Voltage dependence of the nonequilibrium fluctuations of oscillator position (nor- 
malized to those of the uncoupled oscillator) at various temperatures; (a) for oscillator mass 
m = 10~ 12 g, (b) for oscillator mass m = 10 _20 g. 

FIG. 2. Temperature dependence of the nonequilibrium fluctuations of oscillator position (nor- 
malized to those of the uncoupled oscillator) for bias voltage 10 V. 



FIG. 3. (a) Temperature dependence of the nonequilibrium fluctuations of oscillator position 
(normalized to those of the uncoupled oscillator) for bias voltage 10~ 6 U; (b) Magnified low tem- 
perature part of this dependence. 



FIG. 4. (a) Temperature dependence of the nonequilibrium fluctuations of oscillator position 
(normalized to those of the uncoupled oscillator) for bias voltage 10~ 4 V; (b) Magnified low tem- 
perature part of this dependence. 
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